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Challenges to our understanding of radio relics: 
X-ray observations of the Toothbrush cluster 



o 

(N 



o 
u 

6 



> 
CO 

cn 
wo 

cn 
o 
cn 



X 



MA 02138, USA 



UK 



G. A. Ogrean^"*", M. Briiggen^, R. J. van Weeren^, H. Rottgering^ 
J. H. Croston^ M. Hoeft^ 

^Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany 

Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambria _ 
^Leiden Observatory, Leiden University, 2300 RA Leiden, Netherlands 
^School of Physics and Astronomy, University of Southampton, Southampton, S017 IB J, 
^ Thiiringer Landessternwarte, Sternwarte 5, 01118 Tautenburg, Germany 



Accepted xxx xxxx xx. Received xxx xxxx xx; in original form xxx xxxx xx 



ABSTRACT 

The cluster IRXS J0603. 3+4214 is a merging galaxy cluster that hosts three radio 
relics and a giant radio halo. The northern relic, the Toothbrush, is 1.9-Mpc long and 
has an unusual linear morphology. According to simple diffusive shock acceleration 
theory, its radio spectral index indicates a Mach number of 3.3 — 4.6. Here, we present 
results from a deep XMM- Newton observation of the cluster. We observe two distinct 
cluster cores that have survived the merger. The presence of three shocks at or near the 
locations of the radio relics is confirmed by density and temperature discontinuities. 
However, the observation poses several puzzles that challenge our understanding of 
radio relics: (i) at the Toothbrush, the shock Mach number is not larger than 2, in 
apparent conflict with the shock strength predicted from the radio spectrum; (ii) at 
the Toothbrush, the shock front is, in part, spatially offset from the radio emission; (iii) 
at the eastern relic, we detect a temperature jump corresponding to a Mach number 
of approximately 2.5, but there is no associated surface brightness discontinuity. We 
discuss possible explanations for these findings. 

Key words: galaxies: clusters: individual: IRXS J0603. 3+4214 - X-rays: galaxies: 
clusters - shock waves 



1 INTRODUCTION 

In recent years, combined radio, X-ray, and optical observa- 
tions have conclusively shown that giant radio relics - dif- 
fuse, steep-spectrum radio objects observed in the periph- 
eries of s ome clusters - are asso c iated with cluster merg- 
ers (e.^. iGiacintucci et aL 20081: iFinoguenov et aL 1 I2OIOI : 



lMarkevitchll2010l : I van Weeren et al.ll201li v Observations, as 
well as cosmological simulations, suggest that relics trace 
shock waves triggered into the intracluster medium (ICM) 
during major merger events. The basic physical picture is 
that as shock fronts propagate through the ICM, they accel- 
erate particles to relativistic energies and compress magnetic 
fields. The accelerated particles gyrate around the magnetic 
fields lines, consequently emitting synchrotron radiati on ob- 
servable at radio frequencies (e.g., [Fkretti et aI]|2012l V 

The details of acceleration at merger shocks are, how- 
ever, more complex and less well understood. An impor- 
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tant question that still awaits answering is whether shocks 
inject thermal electrons into the cosmic-ray (CR) electron 
population believed to be responsible for the observed syn- 
chrotron emission, or instead they re-accelerate pre-existing 
CR electrons trapped in the ICM magnetic field, which itself 
is frozen-in to the thermal plasma. For strong shocks with 
Mach numbers larger than a few, as found in supernova rem- 
nants, the acceleration process, believed to be diffusive shock 
acceleration (DSA), is very efficient and injection dominates. 
But in the low Mach-number regime [M < 5 for merger 
shocks), the particle acceleration efficiency is too low for 
simple injection to explain the observed radio brightness. At 
such shocks, it appears essential to consider re- acceleration 
of pre-existing CR electrons, as theoreti cal considerations 
and nume rical simulations indicate (e.g. iKang et"aL I I2OI2I : 
IPinzke eTa l. 2013). Yet, observational evidence supporting 
this hypothesis is scarc e at best. 

iKang et al.1 (|2012l ) simulated the particle acceleration 
at the relic in the galaxy cluster CIZA J2242.8+5301. 
They showed that models with pre-existing CRs can ex- 
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plain the relic by a relatively weak shock of Mach num- 
ber 2. However, the injection spectral inde^^ at the relic 
is -0.6 ± 0.05 (jvan Weeren et al.1 12OI0I ). In the DSA test 
particle regime, the injection spectral index relates to 
the shock Mach number via = {2a - 3)/(2a + 1) 

(jBlandford Eichlerl [TqStI ) . Therefore, a Mach number of 
2 would be much lower than the M ~ 4.5 predicted by 
the radio injection spectral index. From the observational 
side, one way to help elucidate the nature of the acceler- 
ating particles is by comparing the radio-predicted Mach 
numbers of putative shocks at radio relics with the ac- 
tual Mach numbers determined from X-ray observations. 
A current observational challenge is posed by the periph- 
eral location of most relics, which makes the shocks dif- 
ficult to detect in X-ray. Consequently, very few shocks 
have been detected so far a t radio relics ('Pinogueno v et all 
2OIOI: [M acario et a l.. l2qil | : [ Ak amatsu & Kawahara I2OIII : 
Qgrean Bruggenll2012l : iBourdin et al.ll2013l ). For the few 
shocks discovered, the X-ray-derived Mach numbers are con- 
sistent with the radio predictions. However, the small sample 
and the large measurement uncertainties make it impossible 
to draw definite conclusions about particle acceleration at 
merger shocks. 

Here, we present results from a deep XMM-Newton 
observation of IRXS J0603. 3+4214, a merging cluster at 
z — 0.225 that hosts thre e radio relics and a giant radio halo 
(Ivan Weeren et aDl2012l ). The northern relic, the largest of 
the three at a length of 1.9 Mpc, has an unusually linear 
morphology, with a broader part to the west; hence its nick- 
name, the Toothbrush. The other two relics have lengths 
of approximately 900 and 200 kpc, and ar e located east and 
southeast of the merger (see, e.g. . fig. 3 of Ivan Weeren et al.l 
I2OI2I V Numerical simulations bv lBriiggen et al. (|2012l l have 
shown that the linear morphology of the northern relic can 
be reproduced by a triple merger with mass ratios 1:1:0.07, 
although the simulated system failed to explain the eastern 
and southeastern relics. Radio observations indicate an in- 
j ection spectral index at the northern relic of —0.6 to —0.7 
van Weeren et al]l2012l ). Therefore, in the DSA test parti- 
cle regime, the Mach number at the Toothbrush is predicted 
to be M = 3.3 — 4.6. Here, we test this prediction by lo- 
cating and characterizing the shocks in the ICM of IRXS 
J0603.3+4214. 

The paper is organized as follows: Section [2] presents the 
observations and the data reduction. Section [3] details our 
analysis, while in Section 2] we analyse the effect of system- 
atic uncertainties on the spectral measurements. Sections [5] 
and [6] discuss the results. A summary of our findings is given 
in Section [71 

We assume a flat ACDM universe with Hq = 70 
kms"^Mpc"\ Qm = 0.3, and Qa = 0.7. At the redshift 
of the cluster, 1 arcmin corresponds to 217 kpc. Through- 
out the paper, quoted errors are la statistical errors, un- 
less stated otherwise. The normalizations of all the spectral 
components are given in default Xspec units (i.e., cm~^ for 
APEC components, and photons keV~^ cm~^ s~^, measured 
at 1 keV, for power-law components) per square arcmin. All 
temperatures are in units of k eV. Metallicities are expre ssed 
using the abundance table of I Anders Grevess i (Il989l l. 



Table 1. Residual soft proton (ReSP) level for the MOS and pn 
event files, and clean exposure times. We list the ReSP calculated 
in the hard band, i.e. 8-12 keV for MOS, and 10 - 14 keV for 
pn. T^sp < 1-2 i ndicates an event file wi th a ReSP level that is 
likely negligible (|Kuntz SnowdenllSoOsI ). 





MOS 1 


MOS 2 


pn 




0.98 ±0.021 


0.98 ±0.019 


1.14 ±0.024 


Exp (ks) 


72.0 


71.7 


61.8 



2 OBSERVATIONS AND DATA REDUCTION 

IRXS J0603.3+4214 - hereafter, the Toothbrush cluster - 
was observed with the XMM-Newton EPIC cameras on Oc- 
tober 3-4, 2011. The observation used the medium filter, and 
had a total exposure time of 82 ks. We analysed the data 
using the XMM-Newton Extended Source Analysis Software 
(esas) integrated in the Scientific Analysis System (sAs) ver- 
sion 12.0.1, and the latest calibration files as of December 
15, 2012. The data reduct ion steps are t he sa me as those 
described in more detail in I Ogrean et al.l (|2012l V 

In summary: Raw MOS and pn event files were first 
filtered of soft proton flares. The filtered event files had ex- 
posure times of 72 (MOS) and 62 ks (pn). We evaluated the 
level of residual soft proton (ReSP) in th e flare- filtered event 
files using the ESAS-tailored approach (iKuntz Snowden 
"2008") of the algorithm presented by iDe Luca Molendil 
(2004). The method is based on quantifying the level of 
ReSP as the ratio of the count rates inside (in-FOV) and 
outside (out-FOV) the field of view (FOV) in the hard en- 
ergy band (8 - 12 keV for MOS, 10 - 14 keV for pn). The 
in-FOV region is chosen to be an annulus at the edge of the 
FOV, where there is typically little ICM emission for clusters 
at z > 0.2. The region selection, along with the rapid drop 
in the detector effective area at high energies, means that 
the in-FOV count rate has es sentially no contribution from 
X-ray photons. According to iKuntz SnowdenI (|2008h , in- 
FOV/out-FOV ratios below 1.2 correspond to event files in 
which the ReSP level is likely negligible. The in-FOV/out- 
FOV ratios (7^sp) for each of the three flare-filtered EPIC 
event files are listed in Table [1] 

The instrumental background was modelled using data 
from the corners (unexposed pixels) of the detectors, and 
filter-wheel closed datasets with hardness ratios and count 
rates similar to those measured during the observation. 
Point sources were detected using ESAS-specific routines, and 
excluded from the analysis. Out-of-time (OoT) events were 
also subtracted from the pn images and spectra. CCD #6 
of MOSl became unoperational after a micrometeorite hit 
in 2005, and is automatically disregarded by ESAS. 



3 DATA ANALYSIS 
3.1 Imaging analysis 

Figure [1] presents the EPIC 0.5 — 4 keV surface brightness 
map of the Toothbrush cluster. The map was instrumental 
background-subtracted, vignetting-corrected, and binned by 
2. The sub-band 1.2 — 1.9 keV, which is dominated by fluo- 
rescent instrumental lines (Al Ka ^ 1.49 keV, Si Ka ^ 1.75 
keV), was not included in the image. The null pixels passing 
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Figure 2. Left: Adaptively-smoothed, vignetting-corrected, and instrumental background-subtracted EPIC surface brightness map, in 
the energy band 0.5 — 4 keV. The adaptive smoothing had a target SNR of 5. Right: Same image, with overlaid WSRT radio contours. 
Contours are drawn at [1,2,4, ...] x 120 /iJy/beam, using a combined radio brightness map created from images at wavelengths 18, 21, 
and 25 cm. 





Figure 1. Vignetting-corrected and instrumental background- 
subtracted XMM-Newton EPIC surface brightness map of the 
cluster, in the energy band 0.5 — 4 keV. The image was binned by 
a factor of 2. The gaps crossing vertically through the centre of 
the image are the result of overlapping MOS and pn CCD gaps. 



Figure 3. Unsharp-masked image emphasizing the surface 
brightness discontinuities. Overlaid are the same WSRT radio 
contours as in Figure O Point sources and zero-exposure pix- 
els are shown in black, to distinguish between shock features (in 
white) and artifacts around no-flux regions. 



in a straight vertical line through the centre of the image are 
caused by the superposition of CCD gaps in the three detec- 
tors. Figure [2] shows the EPIC adaptively smoothed surface 
brightness map in the same energy band. Two distinct clus- 
ter cores are evident in the maps. Both of them have been 
significantly disturbed by the merger. A bridge connects the 
northern, less massive core, to the northwestern part of the 



southern core. The bridge is brighter in the north, closer to 
the smaller core. 

Surface brightness discontinuities are present south of 
the more massive core, and NW of the subcore. They are 
more clearly seen in the unsharp-masked image in Figure O 
The unsharp-masked images were created using two 0.5 — 4 
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Table 2. Best-fit sky background model. Frozen parameters are 
marked with a dagger (f). All parameters are expressed in default 
XSpec units. The normalizations were scaled to correspond to an 
area of 1 sq arcmin. 



Component Tx F A/x 



LHB 


o.ost 


(4.1 ±0.75) X 10-'^ 


GH 


-1 O+0.0032 

"^•-•-^-o.ooso 


(2.2 ±0.086) X 10-6 


HF 


r^ Q7+0.050 

-0.052 


(1.7 ±0.21) X 10-'^ 


CXB 




1.4lt (6.5 ±0.30) X 10-'^ 



keV surface brightness maps smoothed with Gaussians of 
widths 25 and 100 arcsec, by dividing their difference to their 
sum. The large smoothing scale of 100 arcsec removes sub- 
structure from the image, while preserving the global X-ray 
morphology. The smaller smoothing scale of 25 arcsec sur- 
presses Poissonian noise, and reduces small-scale substruc- 
ture and pixel-to-pixel variations (e.g., introduced by the 
GoT events correction, CCD gaps). We are searching for sur- 
face brightness edges on scales of 5 — 10 arcmin (the Tooth- 
brush relic has a length of almost 9 arcmin), much larger 
than the smoothing scales used. We tested the robustness of 
our unsharp-masked images by experimenting with differ- 
ent smoothing; scales an d unsharp-masking procedures (e.g., 
ISanders FabianI l2012l ). The enhanced surface brightness 
edges are visible for various smoothing scales in the range 
of 10 — 200 arcsec, but our choices of 25 and 100 arcsec are 
the best compromise between making the edges clear and re- 
moving substructure. We point out that the same unsharp- 
masking technique has been used bv I Russell et al.l (|2010[ ) to 
highlight substructure and surface brightness discontinuities 
in Abell 2146. The white regions in the map indicate shock 
features, except when they surround excluded point sources 
or zero-exposure pixels (shown in black). The fact that the 
shocks are asymmetric with respect to the merger axis (the 
imaginary line connecting the two cluster cores) indicates a 
merger with non-zero impact parameter. Interestingly, the 
NW shock, which was expected to trace the Toothbrush, is 
visible only along the broadest part of the relic, and extends 
about 700 kpc further to the SW, beyond the radio emission. 

In summary, the surface brightness map shows evidence 
of a merging galaxy cluster involving at least two clusters. 
Two cluster cores have survived the merger, and are con- 
nected by a dense plasma bridge. Shocks triggered during 
the merger are observed as surface brightness discontinu- 
ities N-NW and S-SE of the northern and southern cores, 
respectively. The orientation of the shocks with respect to the 
merger axis suggests a non-zero impact parameter. The N- 
NW shock is visible only near the broadest part of the Tooth- 
brush, and extends more than 0. 5 Mpc beyond it, towards the 
SW. 



3.2 Temperature distribution 

We used the weighted Voronoi tesselations (WVT) binning 
algorithm of Diehl & Statler (2006) to bin the cluster image 
in regions of 3600 counts from source plus sky background, 
starting from the centre of the FGV. We then extracted 
spectra and response files corresponding to each region, and 
binned the spectra to a minimum of 30 counts per bin. Simi- 



Table 3. Best-fit ICM temperatures and normalizations for the 
binned regions in Figure 0] We only list bins for which both 
the temperature and the normalization were constrained. The 
metallicity was frozen to 0.2 solar, using the abundances table of 
I Anders Grevess3 (Il989|) and the photoelectric absorption cross- 
sections of lVerner et al.l |l996h . Values are given in default XsPEC 
units. Normalizations are scaled to an area of 1 sq arcmin. 



Bin Tx A/k 

2 
3 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
21 
22 
23 
24 
25 
26 
27 
28 
29 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 



larly, we also extracted global spectra from a circular region 
encompassing most of the cluster emission. Modelling the 
ICM spectra requires a knowledge of the sky background. 
Therefore, sky background spectra were extracted from a 
partial annulus at the edge of the FGV, where the surface 
brightness is the lowest; this region is least likely to contain 
ICM emission. Like the source spectra, the sky background 
spectra were also binned to a minimum of 30 counts per bin. 
To contrain better the low-energy background components. 



^•U_o.44 ^•-'--0.40 ^ l^ 

8 0+^-^ 2 2+°-°^^ X 10-5 

^•^-0.91 ^-^-0.065 ^ 

5 «+0.63 2 Q+0-12 X 10-5 

^•"-0.50 ^-^-0.10 ^ 

70 + 0.75 -I + 0. 036 y -ipv-4 

'•°-0.73 -^•'-*-0.035 ^ 

7 + 0. 69 -I p^+0.039 -|p)-4 

'•"-0.67 -^-"^-0.038 ^ 

Q 7+0.95 -I + 0. 036 y -ipv-4 

^•'-0.94 -^•^-0.035 ^ 

o 9+1.0 o + 0. 093 y -ipv-4 

"•^-0.75 '-'•'-'-0.091 ^ 

9.1 ± 1.2 3.2+g-g^^ X 10-5 

g + 1.3 2 1 + 0-056 -.Q-4 

y-o_1.2 ^•-'--0.055 X 

'^•3io60 (4.2 ±0.11) X 10-4 

"•^-0.51 ^-"-0.075 X 

7 7+0.56 4 4+0.095 .r.-4. 

' • ' -0.54 ^-^-0.093 X 

9.2+g-|^ (1.9 ± 0.037) X 10-4 

9.6 ±1.1 (4.6 ±0.11) X 10-5 

8-1^0 65 (4.6 ±0.11) X 10-4 

7 8+°-55 1 9+0.038 -I Q-4 
'•^-0.54 -'-•^-0.037 ^ 

8.0lg-^° (5.5 ±0.12) X 10-4 

o n+0."57 o 7+0.080 ^ -|p.-4 

^•^-0.55 '^•'-0.079 X 

8 4+1-3 3 4+0.086 -iQ-4 
^•^-0.66 '^•^-0.093 ^ 

9 2+1-1 g+0.046 -iQ-4 
^•^-0.91 -'-•^-0.045 ^ 

6-Oio 52 (4.8 ±0.15) X 10-5 

o + 1.0 -I c+0.040 -|p,-4 

^•^-0.76 -'-•^-0.034 ^ 

Q -,+1.2 o 0+0. 094 -,p>-4 

^•-'--0.90 '^•^-0.091 ^ 

q -1+0.96 -I p^+0.038 -,^-4 

^•-'--0.95 -'-•^-0.032 ^ 
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8-4io63 (5.0 ±0.12) X 10-4 

gg+0:93 4 7+0.12 X 10-4 

^•"-0.64 ^-'-0.099 ^ 

4 6+°-^^ 6 4+0-46 X 10-6 

^•"-0.73 "-^-0.44 ^ 
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8 3+1-1 . 4+0.14 -IQ-4 

^•^-0.78 ^-^-0.13 ^ 
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^•^-0.60 "-^-0.13 X 1^ 
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^•^-0.70 ^-^-0.15 X 1^ 
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7 f.+0.91 -I 7+0.048 -1^-5 
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Figure 4. Left: Background regions. The XMM-Newton region is shown in red, while the ROSAT regions are shown in blue. Contours 
are X-ray surface brightness contours. Right: Bin number map. Each bin has roughly 3600 ICM plus sky background counts. Gray areas 
correspond to detected point-like sources and to regions outside the central 24 arcmin. 





Figure 5. Temperature map of the cluster. Overlaid are X-ray contours (left) and radio contours (right). Each bin has approximately 
3600 counts from source plus sky-background. 
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we also used three ROSAT Ah-Sky Survey (RASS) spectra 
extracted from regions near the cluster. The sky background 
regions and the binned cluster map are shown in Figure [H 

We fitted all the spectra in parallel using XSpec 
V. 12.7.1. Before fitting, the corresponding instrumental 
background spectra were subtracted from the source and 
sky background spectra. The sky background was mod- 
elled as the sum of Local Hot Bubble (LHB), Galactic 
Halo (GH), Hot ForegrouncO (HF; e.g., ISimionescu et al.l 
I2OIII ). and Cosmic X-ray Background (CXB) emission. 
For the thermal components (LHB, GH, HF), the abun- 
dances were fixed to solar and the redshifts to zero, 
as cu stomary. The LHB temperature was f i xed t o 0.08 
keV (ISidher et al.1 Il996l : iKuntz SnowdenI I2OOOI ). The 
CXB was modelled as a power-law coniponen t, with a 
frozen index of 1.41 (|De Luca Molendil l2004l ). Because 
no point sources were subtracted from the RASS spectra, 
the RASS CXB normalization was fixed to 8.85 x 10~ ^ 
photons keV~^ cm"^ s~^ arcmin"^ (jMoretti et al.1 l2003l ). 
For describing ICM emission, we added to the sky back- 
ground model a single-temperature thermal component with 
a fix ed redshift of 0.225 and a fixed metallicity of 0.2 solar 
(e.g. , iLeccardi Molendilliooi ) . The absorbed components, 
i.e. all except the LHB, were multiplied by a photoelectric 
absorption model with a free X-ray column density; an ex- 
ception are the background spectra, for which the column 
density was fixed to 2.14 x 10^^ cm~^ - the Galactic Hi 
column density listed in the L eiden/Argentine/Bonn (LAB) 
Survey (|Kalberla et al.l [20051 ) in a 0.5- degree circle around 
the cluster - because of degeneracy between the background 
model normalizations and the absorption. The fluorescent 
instrumental lines of Al Ka (E ^ 1.49 keV) and Si Ka 
(E ^ 1.74 keV) vary in intensity with time and across each 
of the EPIC detectors, and therefore they are not included 
in the instrumental background model. We modelled them 
separately using two zero- width Gaussian components with 
free normalizations, but with central energies coupled be- 
tween spectra. 

The best-fit background model is summarized in Ta- 
ble [21 Table [3] shows the ICM temperatures and normal- 
izations for all binned regions for which we were able to 
constrain the spectral parameters. The fit had a x^/d.o.f. = 
12380.43/12110. Figure [5] presents the resulting temperature 
map, with overlaid X-ray and radio contours. 

The duster has an average temperature of 7.8 dz 0.13 
keV, and Lx, 0.5-7 keV = 1-9 x lO'^^ erg s~^ (Lx,o.i-2AkeV = 
7.7 X 10^^ erg s~'^ ). Along the Toothbrush, the temperature 
decreases from approximately 8 keV to the SW, to about 
2 keV to the NE. E-SE of the merger, there is a region of 
hot plasma with temperatures of around 11 keV. For the bins 
shown in Figure\^ the average temperature is 8.1 keV, with a 
standard deviation of 1.6 keV. Most temperatures are within 
la of the average when statistical errors are considered. 



4 SYSTEMATIC UNCERTAINTIES 

There are three sources of systematic errors on the measured 
ICM parameters: the X-ray column density in the sky back- 
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^ The Toothbrush cluster is located at a Galactic latitude of ? 
degrees, where hot foreground emission is possibly present. 
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Figure 7. Temperature map, with statistical and systematic er- 
rors taken into account. All bins with temperatures within Icr 
of the average are shown in orange. Overlaid are the same radio 
contours as in Figure [2] 



ground region, systematic uncertainties on the sky back- 
ground parameters, and uncertainties on the instrumental 
background model. 

In a circular region of radius 1-degree around the clus- 
ter, there are seven Galactic Hi measurements, with corre- 
sponding Galactic hydrogen column densities in the range 
(1.86 — 2.47) X 10^^ cm~^. The weighted average column 
density is 2.15 x 10^^ cm~^, essentially equal to the column 
density measurement at the position closest to the cluster. 
To analyze the effect of an uncertain X-ray column den- 
sity in the sky background region, we refitted all the spec- 
tra, once with the sky background column density fixed to 
1.86 X 10^^ cm~^ (the lowest A^h in the region), and once with 
it fixed to 2.47 x 10^^ cm"^ (the highest A^h value). Using 
each of the two new sets of best-fit background parameters 
and the corresponding statistical errors, we refitted the ICM 
spectra while varying (one by one) the sky background pa- 
rameters by their la statistical errors; the sky background 
model was fixed in these fits. The resulting minimum and 
maximum ICM parameters for each bin define the ranges for 
these parameters when systematic uncertainties on the sky 
background parameters and column density are taken into 
account. 

Figure [6] shows the best-fit temperatures listed in Ta- 
ble O with overlaid la statistical error bars, and the ranges 
of systematic uncertainties. The statistical errors dominate, 
so A^H uncertainties and la variations on any of the back- 
ground parameters have only a small effect on the best-fit 
ICM temperatures. The same is also true for normalizations 
(not shown in Figure[6]), for which the effect is even stronger. 
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Figure 6. Top: Statistical and systematic errors on the best-fit temperature measurements. Charcoal bars show the best-temperatures 
listed in Table (3] while bisque bars show the minimum and maximum temperatures when systematic uncertainties on the sky background 
parameters are taken into account. Bottom: ICM temperature uncertainties introduced by a 5% variation in the normalization of the 
instrumental background. Charcoal bars show the best-temperatures listed in Table[3] while bisque bars show the minimum and maximum 
temperatures obtained when the instrumental background level is varied by ±5%. In both plots, error bars are la statistical errors. 



We evaluated separately the effect of a possible 5% error 
on the normalization of the quiescent particle background. 
For each bin spectrum (see Figure |4]), the normalizations 
of the MOS and pn spectra were simultaneously increased 
and then decreased by 5%. The sky background was fixed to 
the model summarized in Table O and the spectra were fit- 
ted again to obtain new best-fit ICM temperature and nor- 
malization, and X-ray column density. The minimum and 
maximum ICM parameters yielded by these fits are shown 
graphically in Figure [6] For most bins, the temperature un- 
certainty introduced by a potential 5% error on the QPB 
normalization is smaller than the statistical error. Excep- 
tions are the outer bins (bins 2, 33, 45, and 46), for which 
5% changes in the instrumental background level resulted in 
possible temperature ranges significantly larger than 2a. 

In Figure [7] we show again the temperature map, includ- 
ing the uncertainties introduced by statistical and system- 
atic errors. Bins with temperature ranges consistent, within 
1(7, with the average temperature of 8.1 keV are all shown in 
orange. Bins 2, 5, and 41, which have temperatures markedly 
different from the average, are shown in purple, green, and 
dark red, respectively. 



5 SHOCKS 

To examine the surface brightness discontinuities identified 
in the unsharp-masked image, we extracted EPIC instru- 
mental background-subtracted surface brightness profiles in 
the sectors shown in Figure (8] The annuli in each sector 




Figure 10. The adaptively smoothed image shown in Figure [2] 
with overlaid regions showing the sectors used for extracting the 
surface brightness profiles in Figures fTTI and fT2] Dotted lines show 
the best-fit positions of the density discontinuities. Overlaid are 
the same radio contours as in Figure [21 
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Figure 8. The same unsharp-masked and adaptively smoothed images as shown in Figures [2] and [S] with overlaid regions showing the 
sectors used for fitting the surface brightness discontinuities. Dotted fines show the best-fit positions of the density discontinuities. The 
surface brightness profiles and best-fit models are shown in Figure [9] 
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Figure 13. Regions used for calculating the post- and pre-shock 
temperatures at the N (red), S (green), SE (orange), and E (pur- 
ple) shocks. Radio contours are the same as in Figure [21 



were chosen to have a minimum signal-to-noise radio (SNR) 
of 10. The resulting profiles are shown in Figure [9l Each 
profile was fitted with a broken power-law density model 
in a region around the putative shock front radius, plus a 
constant describing the sky background. The density models 



can be summarized as: 

/ ^ \ «2 

722 = C no I I , for r ^ rshock 

\ ^shock / 

f r 

Til — TiQ \ j , otherwise. 

V '^shock / 

In the equations above, n is the electron number density, C 
is the shock compression, r is the distance from the centre of 
the sector, rshock is the shock radius, and a is the power-law 
index. Subscripts 1 and 2 stand for pre-shock and post-shock 
parameters, respectively. The emissivity, proportional to the 
density square, is integrated along the line-of-sight, assum- 
ing spherical symmetry, to calculate the surface brightness 
model. The sky background value was calculated from the 
region shown in Figure [H and kep t fixed in the fit. F itting 
was done using PROFFIT vl.]0 (lEckert et al.ll201lh . The 
best-fit models are shown in Figure [9l 

As expected from the unsharp-masked image, we find 
surface brightness discontinuities in each of the two sectors. 
To the north, the assumed underlying density model indi- 
cates shock compression by a factor of 2.0lg:l^; the fit had 
X^/d.o.f. = 66.2/57. To the south, the compression factor 
is l.Tto n, with xV^-oi. = 85.5/57. The Mach numbers 
at the northern and southerns discontinuities are therefore 
relatively low, 1.71q;42 l-^^o.o86 5 respectively. 

Spherical symmetry is one of the main assumptions be- 
hind our surface brightness models. However, it is a rather 
crude approximation, given the clearly disturbed cluster 
morphology. We can minimize the effects introduced by de- 
viations from spherical symmetry by reducing the opening 
angle of each sector. Hence, we divided each of the sectors 
in Figure [8] into two, and searched for surface brightness 

^ http://www.iasf-milano.inaf.it/ eckert/newsite/Proffit.html 
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Figure 9. Top: Instrumental-background-subtracted, exposure-corrected surface brightness profiles across the northern (left) and south- 
ern (right) sectors shown in Figure [8] The purple lines show the instrumental background in the two sectors. Each bin has a SNR of 
at least 10. Bottom: Surface brightness profiles around the density discontinuities, with best-fit models. The models assume sphericity 
and an underlying density profile described by two power-laws with a jump at the shock front. The sky background value was fixed 
to 4.4 X 10~^ counts s"-*^ arcmin"^, the average background value calculated in the XMM-Newton region shown in Figure 0] For the 
northern profile, black vertical lines (dotted) mark the inner and outer boundaries of the Toothbrush relic, and the green line (dot-dash) 
marks the outer boundary of the N-NW radio extension. For the southern profile, black vertical lines (dotted) mark the inner and outer 
boundaries of the small (~ 200 kpc) SE relic, while the green line (dot-dash) marks the southern boundary of the radio halo. In both of 
the bottom plots, the best-fit shock radii are marked with red lines (dashed). The surface brightness jumps correspond to Mach numbers 
of 1.7l°;^^ (north) and 1.51°;°^^ (south). 



jumps in the resulting profiles. The new sectors and profiles 
are shown in Figures [T0lll2l 

To the north, across the broadest part of the Tooth- 
brush relic, a very weak density discontinuity is present 
at the outer boundary of the small N-NW radio exten- 
sion. This discontinuity corresponds to a Mach number of 
1.3 ±0.45, therefore consistent with 1. To the NW, the Mach 
number inferred from the best-fit density jump is higher, 
l-9io.42 5 although the errors are also large. Both the N and 
NW fits were satisfactory, with x^/d.o.f. = 42.9/35 and 
X^/d.o.f. = 37.9/27, respectively. 

To the south, we divided the larger sector in Figure [8] 
such that one profile crosses the SE relic, while one crosses 
the edge of the radio halo. The SE profile is excellently fitted 
by our model, with x^/d.o.f. = 23.8/23 and a best-fit den- 
sity jump of l-7lo'i8 5 corresponding to a Mach number of 
1-^-0. 13- Modelling the S profile indicates a Mach number 
of 1.3lQ jg7, although the fit is much poorer (x^/d.o.f. = 
70.7/34)'. 

We further examined the detected surface brightness 
discontinuities by extracting spectra in regions tangential 
to the best-fit shock front radii. These regions are shown in 
Figure [13] To the N-NW, the difference between the best- 
fit shock radii in the two sectors shown in Figure [8] is only 



6.1 arcsec, essentially the spatial resolution limit of XMM- 
Newton. In order to improve the SNR of the spectra ex- 
tracted N-NW of the Toothbrush, we have selected two par- 
tial annuli that cover the full width of the N-NW sector in 
Figure [H] The boundary between these partial annuli, i.e. 
the shock radius, is chosen to be the average of the best- 
fit radii for the narrower N and NW sectors, rshock,N = 5.0 
arcmin. The pre-shock and post-shock spectra were fitted 
with the same model as the bin spectra, with the sky back- 
ground parameters fixed to the best-fit values in Table (2] 
The X-ray column density was fixed to 1.97 x 10^^ cm~^, 
the average value of bins 3 (A^h = 1-98 x 10^^ cm~^) and 
14 (A^H = 1.97 X 10^^ cm"^), which overlap with the se- 
lected post-shock and pre-shock regions. The best-fit post- 
shock and pre-shock temperatures are 7.01q'4q and 7.1^2 4 
keV. We performed the same systematic error analysis as 
described in the previous section. The temperature ranges 
resulting from la changes in the sky background and ±5% 
changes in the QPB are 6.6 — 7.8 and 6.5 — 7.5 keV, respec- 
tively, for the post-shock region, while for the pre-shock re- 
gion they are 4.3 — 21.8 and 4.1 — 11.1 keV, respectively. The 
post-shock temperature is very well- constrained. Although 
the uncertainties on the pre-shock temperature are signifi- 
cantly larger, they allow us to set a lower limit of approxi- 
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Figure 11. Top: Instrumental-background-subtracted surface brightness profiles across the NW (left) and N (right) sectors shown in 
Figure [To] The purple lines show the instrumental background in the two sectors. Each bin has a SNR of at least 10. Bottom: Surface 
brightness profiles around the density discontinuities, with best-fit models. The models are the same as those used for the profiles shown 
in Figure O For the N profile, black vertical lines (dotted) mark the inner and outer boundaries of the Toothbrush relic, and the green 
line (dot-dash) marks the outer boundary of the N-NW radio extension. In both of the bottom plots, the best-fit shock radii are marked 
with red lines (dashed). The surface brightness jumps correspond to Mach numbers of I.Q^q'^ (NW) and 1.3 ± 0.45 (N). 



mately 4 keV. This implies a Mach number M < 1.9, con- 
sistent with the value derived from the surface brightness 
jump. 

As an exercise in caution, we ignored the shock front lo- 
cation determined above, and extracted an additional spec- 
trum from a "pre-shock" region tangent to the outer edge 
of the broadest part of the Toothbrush (Figure [15)), rather 
than to the N-NW radio extension; naively, the outer edge 
of the relic is the expected shock front location. However, 
the best-fit temperature in this region was similar to the ac- 
tual pre-shock temperature, 5.3l};^ keV. Even when possible 
systematic errors were taken into account, the temperature 
was still above 3 keV, ranging between 3.6 and 9.4 keV. 
Hence, the Mach number at the Toothbrush relic, inferred 
both from the surface brightness jump and from the tem- 
peratures on both sides of the shock front, is not larger than 
approximately 2. 

To the S and SE, we extracted spectra from partial 
annuli of width 1.5 — 2 arcmin, tangential to the shock 
front positions shown in Figure [TOl The best-fit post-shock 
and pre-shock temperatures were 9.9 ±0.92 and 8.ll2'o 
to the S, and lO.OtiJ and QAt\l keV to the SE. In the 
S regions, the poor count statistics do not allow us to de- 
tect a possible temperature discontinuity. In the SE regions, 
the temperature jump corresponds to a Mach number of 
1.61q'33, consistent with the Mach number determined from 
the surface brightness jump. When systematic errors were 
taken into account for the two SE regions, the best-fit tem- 



peratures varied between 9.5 — 10.5 keV for the post-shock 
region, and 5.3 — 7.9 keV for the pre-shock region; therefore, 
the temperature jump is certain, and the Mach number mea- 
surement across the SE relic is a robust result. 

Across the eastern relic, we searched for shock evidence 
by measuring the temperature in two partial annuli tangent 
to the outer edge of the relic. The temperature decreases 
from S.7t\l keV in the more central region, to 3.31q'57 keV 
in the outer region. This jump indicates a shock front of 
Mach number 2.4lQ'3g. Both the post-shock and pre-shock 
temperatures are well-constrained, even when systematic er- 
rors are taken into account; they range between 7.7 — 9.8 
and 2.7 — 4.1 keV, respectively. However, as can be seen in 
Figure [M] there is no surface brightness discontinuity any- 
where near the relic. A /3-model is an acceptable fit to the 
profile, with x^/d.o.f. = 35.0/25. The data is fitted better 
by a broken power-law model, x^/d.o.f. = 25.2/23, yet al- 
though there is a kink in the profile (the power indices are 
Oil — l-S^o 047 ^2 = 3.2I0 35)' "there is no actual jump 
(the shock compression is C = l-O^o.ii)- 

In summary: Surface brightness and temperature dis- 
continuities indicate the presence of three, possibly four, 
shocks fronts in the ICM, towards N-NW, SE, and E (and 
perhaps towards S). The N-NW shock is about 1 arcmin 
ahead of the Toothbrush's outer edge, extends more than 
half a Mpc beyond the relic 's W tip, and has a Mach num- 
ber significantly below that predicted by the radio spectral 
index. At the eastern relic, we detect a temperature jump 
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Figure 12. Top: Instrumental-background-subtracted surface brightness profiles across the SE (left) and S (right) sectors shown in 
Figure [To] The purple lines show the instrumental background in the two sectors. Each bin has a SNR of at least 10. Bottom: Surface 
brightness profiles around the density discontinuities, with best-fit models. The models are the same as those used for the profiles shown 
in Figure [9] For the SE profile, black vertical lines (dotted) mark the inner and outer boundaries of the SE relic. In both of the bottom 
plots, the best-fit shock radii are marked with red lines (dashed). The surface brightness jumps correspond to Mach numbers of l.S^Q^Jg 
(SE) and 1.3lg:J°7 (S). 



that corresponds to a Mach number of approximately 2.5, 
hut there is no associated surface brightness discontinuity. 
The surface brightness jump could be masked by projection 
effects, as was also suggested fo r the shock front in Abell 665 
iMarkevitch & Vikhlinir\\200i ). The int egrated spectral in- 
dex o f the eastern relic is a = — l.OdbO.2 Hvan Weeren et all 
inconsistent with the X-ray- derived Mach number if 
one assumes an injection spectral index ainj = a + 0.5; since 
the A^(ainj) function diverges for ainj = —0.5. The SE relic 
has a radio -predicted Mach number of 4.6 di 1.1^ also signif- 
icantly above the X-ray result. 



6 THE TOOTHBRUSH PUZZLE 

The results presented above pose several challenges to our 
understanding of radio relics: 

• Not all relics are associated with surface brightness dis- 
continuities. 

• The shock front near the Toothbrush is spatially offset 
from the relic, both in distance and position angle from the 
cluster centre. 

• The Mach number of the shock is significantly lower 
than the Mach number predicted by the injection spectral 
index. 



1.0 - 




-1.0 -0.5 0.0 0.5 1.0 

Mpc 



Figure 15. Additional region used for measuring the temperature 
ahead of the outer relic edge. This region would be the expected 
pre-shock region. The best-fit temperature is 5.3^^"^. Contours 
are the same radio contours as in Figure O 



Below we discuss each of these issues. 
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Figure 14. Instrumental-background-subtracted surface brightness profile in a sector across the eastern relic. The centre of the sector 
is the same as the centre of the S and SE sectors. Each bin has a SNR of at least 15. Blue lines show the best-fit /3-model (left; 
X^/d.o.f. = 1.4) and broken-power law model (right; x^/d.o.f. = 1.1). Red dashed lines mark the position of the outer edge of the relic; 
the same putative shock front radius was used to separate the partial annuli used for measuring the temperature jump across the relic. 



6.1 No surface brightness jump 

Although most large-scale shocks are associated with sur- 
face brightnes s discontinuities, there are als o exceptions. 
In Abell 665, iMarkevitch k VikhlininI (l20Qlh identified a 
shock front in the temperature distribution map, but no 
corre sponding surface brightnes s dis continuity. At the Coma 
relic JOgrean Briiggenl ( 2012h and lAkamatsu et al ] (I2OI3I ) 
found a shock of Mach number 2 based on the temperature 
jump, but no clear surface brig htness discontinuity is see n 
in the surface brightness profile (jOgrean & Briiggen 2012"). 

At the eastern relic in the Toothbrush cluster, the tem- 
perature jump across the relic's outer edge implies a Mach 
number of 2A'tQ]s6- Using the Rankine-Hugoniot jump con- 
ditions, the Mach number corresponds to a density jump by 
a factor larger than 2. However, this jump is not reflected 
in the surface brightness profile; spherical symmetry implies 
a shock compression consistent with 1, with a possible kink 
in the underlying density profile. The assumption of spheri- 
cal symmetry in such a disturbed system is clearly not valid. 
Therefore, the lack of a surface brightness jump is not equiv- 
alent to the lack of a density jump. Even if a density jump 
was present across the shock, it could potentially be masked 
by deviations from sphericity and substructure along the 
line of sight. Their effect on the temperature measurements 
would be less pronounced. 

6.2 Relic-shock offset 

As also mentioned in Section [1] radio relics have only been 
discovered in merging galaxy clusters, pointing to an asso- 
ciation between the two. It is believed that shocks triggered 
during cluster mergers can accelerate electrons to relativis- 
tic energies, and these accelerated electrons gyrate around 
magnetic field lines, emitting synchrotron radiation at ra- 
dio frequencies. The acceleration mechanism is believed to 
be diffusive shock acceleration. Assuming direct acceleration 
at the shock, the shock front location is expected to be at 
the outer edge of the relic. As shown in the previous sec- 
tion, this is not the case at the Toothbrush. The shock front 
position does not coincide with the outer edge of the Tooth- 
brush, with a confidence of approximately 7a. Even more 



puzzling, while the relic extends towards the NE, the shock 
front extends towards the SW, although the two appear to 
be roughly parallel. There are several possible explanations 
for the offset: 

• projection effects; 

• changes in magnetic field between the relic and the 
shock front; 

• a pre-existing population of CR electrons; 

• a non-detected density discontinuity at the outer edge 
of the Toothbrush. 

Projection effects are unlikely though to explain the 1- 
arcmin offset between the shock and the relic. If relics are 
created by direct acceleration at the shock, the shock is co- 
spatial with the outer edge of the relic. Assuming the shock 
and the relic are spherical caps, it would be possible to de- 
tect the shock front on the projected surface of a broad 
relic. However, it would be impossible to project the de- 
tected shock front ahead of the relic. 

From equipartition argume nts, the magnetic field at the 
Toothbrush relic is 7 - 9 /xG (jvan Weeren et al.1 lioTi ). If 
the magnetic field strength dropped between the observed 
locations of the shock front and the Toothbrush, then syn- 
chrotron emission might not be observable in that region, 
creating an offset between the shock and the relic. However, 
magnetic fields are compressed at the shock (and possibly 
also amplified by downstr eam turbulence, a s seen at higher 
Mach number shocks; e.g. Jinoue et al.ll2009h , so it is unclear 
how the magnetic field could drop suddenly in front of the 
Toothbrush. 

Another possibility is that the relic is the result of re- 
acceleration of a pre-existing CR electron population. One 
problem with the simplest DSA model is that the parti- 
cle acceleration efficiency is too low to explain t he radio 
brightness of relics (e.g., iKang et al.ll20Q7l . I2QI2I ). Conse- 
quently, the model has been challenged by numerical simu- 
lations showing that a pre-existing population of CR elec- 
trons makes it easier to naturally exp l ain the observed ra - 
dio properties (e.g., iKang et al.ll2012l : IPinzke et al.ll2Q13l ). 
CRs are expected to be formed by, e.g., la rge scale structure 
formation shocks (e.g., IRvu et aU l2003h , turbulence (e.g.. 
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iBmnetti et"al]l200lh . hadronic collis ions of CR protons wit h 
thermal protons in the IC M (e.g., IPolag k EnfilinI l2Q00h . 
active galactic nuclei (e.g., iGeorge et al.ll2008l ). and super- 
novae (e.g., Blasi 2011), and they will remain trapped in the 
cluster's potential well. If a pre-existing CR population was 
present only at the location where the relic is seen, then as 
the shock observed to the N-NW crossed this region, it would 
have re-accelerated the CRs to energies where they become 
visible in the radio. On the other hand, beyond the location 
of the CR population, the efficiency of particle acceleration 
from the thermal pool into the non-thermal CR population 
would be low, so no radio emission would be detected. 

An unusual situation would be the presence of a density 
edge at the relic, in addition to the shock detected further 
out from the cluster centre. This edge could either be a sep- 
arate shock, or we could be looking at a complex shock sur- 
face for which a surface brightness discontinuity is seen only 
offset from the relic's outer edge. Figures [9] and [11] do not 
show a surface brightness discontinuity at the relic's outer 
edge. Yet, as we discuss in Section [6Jl multiple effects and 
assumptions can mask a potential density discontinuity. 

6.3 Mach number discrepancy 

In the linear test-particle regime, DSA predicts that the in- 
jection spectral index at the shock front, a (J^ oc zy"), is 
related to the shock Mach number, A^, via 



(2a-3)/(2a + l) 



(1) 



for plane-parallel shocks. If particles are accelerated directly 
from the thermal pool and the assumptions of the test- 
particle regime hold, then the radio-derived Mach num- 
ber should be consistent with the Mach number derived 
from X-ray data using the Rankine-Hugoniot jump con- 
ditions. Indeed, this is the case for all the other merger 
shocks confirmed at radio relics. In Abell 3667, the spec- 
tral index near t he outer edge of the radio relic is ap- 
proximately —0.7 (jRottgering et al.lll997l V corresponding to 
a Mach number of approximately 3.3; the X-ray-derived 
Mach number at the relic was found to be 2.4 it 0.77 
(|Finoguenov et al]|2010h from the temperature jump across 
the shock front (we cite the Mach number indicated by the 
temperature jump, as temperature is less affected by pro- 
jection effects than density). In Abell 754, the integrated 
radio spectral index is between —2.0 (for the frequency 
range 0.3 — 1.4 GHz) and —1.8 (for the frequency range 
0.3 — 0.7 GHz) ; if the injection spectral index is ainj = a+0.5 



[zj ; 

(|Miniatill2002l ). the spectral indices at these frequencies cor- 
respond to a Mach number of 1.7 — 1.9, consistent with the 
Mach number of 2.lln'^ d etermined from the temperature 



jump (|Macario et al.ll201lh . In CIZA J2242.8+5301, the so- 
called "Sausa ge" relic has a radio-pr edicted Mach number 
of 4.6 =b 1.1 ( van Weeren et al"]|2010l l. consistent with the 
Mach number of 3.2 di 0.52 determined from Suzaku obser- 
vations near the relic ([Akamatsu Kawaharal [201lh . The 
western relic in Abell 3376 has a spectr al index that cor - 
responds to a Mach number of 2.2 db 0.4 (>Kale et al.ll2012l \ 
which is consistent with the X-ray-der ived Mach number of 
2.9=b0.6 (| Akamatsu k Kawaharall201lh . In Abeh 521, the ra- 
dio spectral index predicts a Mach number of approximately 
2.3 ([Giacintucci et al1l2008l V consisten t with the Mach num - 
ber of 3.4l^ g determined from X-ray (jBourdin et al.|[2013l V 



The injection spectral index at the Toothbrush relic is 
between —0.7 and —0.6, corresponding toA^ = 3.3—4.6. The 
Mach number derived both from the density and the temper- 
ature discontinuities across the shock front is no larger than 
2, significantly lower than the simplest DSA predictions. A 
possible explanation is that any shoc k front is comprised o f 
a distribution of Mach numbers (e.g. JSkillman et al]|2013l V 
Sync hrotron emission is mor e sensitive to high Mach num- 
bers (jHoeft Bruggenll2007l ) , so shock regions of high Mach 
number can introduce an upward bias in the overall Mach 
number derived from radio o bserva tions. 

Moreover, iKang et al.1 (|2012l ) have shown that the 
Sausage re lic, which has a radio-pr edicted Mach number of 
4.6 ±1.1 (|van Weeren et all I2OIOI I. can be also explained 
by a shock front with Mach number of about 2, if a pop- 
ulation of pre-existing CR electrons is introduced in the 
model. While such a low Mach number at the Sausage 
relic ha s been shown not to be co nsistent with obser- 
vations (lAkamatsu Kawaharal l201lh , the conclusions of 
iKang et al ] (I2OI2I ) suggest that one possible explanation for 
the Mach number discrepancy at the Toothbrush relic is 
that the shock re-accelerated a pre-existing population of 
CR electrons. 

Alternative explanations also exist. Projection effects, 
for example, could smooth out both the density and the 
temperature discontinuities. Furthermore, the Mach num- 
ber range of 3.3 — 4.6 was derived from the radio spectral in- 
dex assuming a plane-parallel shock. For oblique shocks, the 
simple equation relating the Mach number and the spectral 
index (Eq.[T]) over-pr edicts the strength of the shock (e.g., 
iKirk k Heavenslll989l ). 

The shock front detected at the Toothbrush relic chal- 
lenges our current understanding of radio relics. It is spa- 
tially offset from the radio relic, and its Mach number is sig- 
nificantly below that predicted from the spectral index at the 
front of the relic assuming DSA in the test particle regime. It 
appears challenging to explain the observed shock character- 
istics using the simplest DSA model, in which CR electrons 
at the relic are accelerated directly from the thermal pool. 
However, the observations alone cannot exclude alternative 
explanations such as projection effects or complex magnetic 
fields, and numerical simulations would be required to test 
the different hypothesis. 



7 SUMMARY 

The Toothbrush cluster {z = 0.225) hosts three radio relics 
and a faint radio halo. The northern relic, the Toothbrush, is 
the most spectacular of the three relics due to its unusual lin- 
ear shape, which might have been caused by a triple merger 
event. The spectral index at the front of the relic suggests 
a Mach number between 3.3 and 4.6. Here, we analyzed a 
deep XMM-Newton observation of the cluster, identified the 
shock fronts present in the ICM, and compared their prop- 
erties with those predicted by radio observations. Below is 
a summary of our main results: 

• We find clear evidence both from imaging and spectral 
results that the Toothbrush cluster is a merging system. The 
cores of the merging clusters survived the merger, and are 
now observed at a separation of 650 kpc, connected by a 
bridge of dense plasma. 
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• The temperature decreases from W to E along the 
Toothbrush. 

• Surface brightness discontinuities are present N-NW of 
the subcluster's core, and S-SE of the more massive core. 

• We find evidence for three, possibly four, weak shock 
fronts. They are found towards N, E, and SE (and possibly 
S), near the positions of previously discovered radio relics, 
and have Mach numbers lower than about 2.5. All three 
Mach numbers are smaller than the Mach numbers predicted 
from the radio spectral indices under the assumptions of di- 
rect CR injection and DSA in the linear test particle regime. 
The shocks in the Toothbrush cluster are the only ones with 
such clear differences between the radio-predicted and X- 
ray-derived Mach numbers, and a larger relic sample with 
very high quality radio and X-ray data would be required 
for a statistical analysis of such occurences. With the advent 
of new radio arrays such as LOFAR and SKA, more relics 
with very high Mach numbers predicted by their spectral in- 
dices could prove instead to be associated with much weaker 
shocks. 

• The N-NW shock, which was expected to trace the 
outer edge of the Toothbrush relic, is unlike any of the other 
shocks discovered so far at radio relics. It is offset by about 
1 arcmin NW from the edge of the Toothbrush, and extends 
by about 700 kpc towards SW, beyond the W tip of the 
relic. Numerical simulations will be required to explain the 
shock's unusual properties in relation to the relic. 
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